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(57) Abstract: The invention relates to a method of reconstructing a radiographic image by combining elemental images. In a 
2£ radiography method whereby an object (2) is examined in numerous steps using a radiation (4) that occupies different positions (4i, 
4j), as does the associated network of detectors (3i, 3j), the elemental images are combined into a complete image without junction 
defects by (i) discretising the object (2) into volumes and (ii) calculating the attenuation in each volume (voxel) (8), in order to obtain 
images of the object at different reconstruction heights. Subsequently, said images are combined together to produce a more precise 
O final image. 
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En ce qui concerns les codes a deux lettres et autre s abr eola- 
tions, se referer aux "Notes explicatives relatives aux codes et 
abreviations "figurant au debut de chaque numero ordinaire de 
la Gazette du PCT. 



(57) Abrege : Dans un proc£d£ de radiographic ou un objet (2) est etudie en plusieurs Stapes par un rayonnement (4) occupant 
diverses positions (4i, 4j) de meme que le rSseau de dgtecteurs (3i, 3j) associe\ la combinaison des vignettes elementaires en une 
image d'ensemble sans dSfaut de jonction est resolue par une discretisation de Fobjet (2) en volumes et un calcul de l'att£nuation 
dans chaque volume (voxel) (8) pour obtenir des images de l'objet a differentes hauteurs de reconstruction ; puis ces images sont 
combiners entre elles pour obtenir une image finale plus exacte. 



WO 2004/059578 



PCT/FR2003/050195 



1 



PROCEDE DE RECONSTRUCTION D'UNE IMAGE RADIOQRAPHIQUE 
PAR COMB INAI S ON DE VIGNETTES SE RECOUVRANT 

5 DESCRIPTION 

L 1 invention concerne un procede de 
reconstruction d'une image radiographique par 
combinaison d'une collection de vignettes se 
10 recouvrant. 

Une application importante, mais qui n'est 
pas exclusive, de 1' invention est 1 'ostfiodensitometrie, 
c'est-^-dire la mesure de density min^rale osseuse 
(DMO) du corps, o\l il est aussi entrepris des mesures 

15 de la composition des tissus nous en distinguant les 
tissus maigres des tissus gras. L'examen peut porter 
sur des regions ^tendues du corps . La DMO s 1 exprime 
comme une masse par unit£ de surface qui correspond & 
la projection, d'aprds des lignes parall§les, de la 

20 masse osseuse sur un plan rapportSe a une surface 
unitaire. La multiplication de la DMO par la surface 
des os donne le contenu mineral osseux ou CMO. Un 
int£r§t de I 1 invention dans cette application sera non 
seulemeht de donner des images meilleures comme on 

25 l'indiquera plus loin, mais des mesures particuli^res 
plus exactes. 

Les images de grande taille en radiographic 
sont f r£quemment obtenues par morceaux, au moyen- de 
vignettes de projection qui" sont prises successivement 

30. en. d£pla<?ant le rayonnement traversant l'objet k 
diffSrentes positions, ainsi que le riseau 
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monodimensionnel ou bidimensionnel de d^tecteurs 
prenant les mesures. L' assemblage des vignettes doime 
ensuite. 1' image recherchSe. 

Le proc£d£ se coraplique dans le cas usuel 
5 d'un rayonneraent divergeant d'un foyer vers le rlseau 
de d£tecteurs, soit en c6ne, soit en un ensemble 
d'6ventails plans et parall£les. La figure 1 montre la 
configuration normale des mesures : le rayonnement 
comprend une source 1 (ponctuelle ou linSaire) qu f on 

10 d^place 1 chaque mesure le long de l'objet 2 ainsi que 
le r£seau de dStecteurs 3 • Les positions prises sont 
not^es la, lb, lc, Id, et 3a, 3b, 3c et 3d. Pour que 
1 ' attenuation du rayonnement soit mesuree §l tout 
endroit de l'objet 2, le faisceau 4 du rayonnement doit 

15 comprendre des portions de recouvrement , assez grandes 
pour que tant point de l'objet 2 soit vu au moins- une 
fois, dans les positions 4a, 4b, 4c et 4d qu'on lui 
fait prendre, et les vignettes de projection du 
rayonnement, dont les positions coincident avec celles 

20 3a, 3b, 3c et 3d que le rlseau de detecteurs 3 prend 
success ivement, pr6sentent de mSme des portions de 
recouvrement. II est done impossible de juxtaposer 
simplement les vignettes pour obtenir 1 1 image globale 
de l'objet, mais on doit au contraire determiner les 

25 positions des portions de recouvrement sur les 
vignettes et faire une synthase du contenu de ces 
portions de recouvrement pour reconstruire 1 ' image . 

Un autre probl£me qui apparalt est celui du 
grossissement des details selon leur distance & la 

30. source 1.. La largeur .de projection de details 5 de 
l'objet .2 . sur le reseau de d£t«cteurs 3 sera 
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proportionnell'ement plus large si les details 5 sont 
plus proches de la source 1- Un rayonnement divergent 
ne permet ainsi ni de juxtaposer facilement des 
vignettes, ni de respecter 1'echelle des details S 
5 l'int^rieur de chaque vignette. 

La figure 2 permet d'expliciter ces 
problSmes. Deux details 5a et 5b espaces verticalement 
se trouvent dans la portion de recouvrement des 
vignettes prises par le r^seau de d^tecteur 3 aux 

10 positions 3a et 3b. Les rayons passant par le detail 5a 
sont distants de l"6cart 6 sur le plan du r£seau de 
d^tecteurs 3, et ceux qui traversent le detail 5b sont 
distants de 1 • 6cart 7 sur le mSme plan ; les hearts 
sont diffSrents entre eux, et diff^rents du d^placeraent 

15 qu'il a fallu imposer au reseau de d^tecteurs 3 entre 
les positions 3a et 3b oil les vues ont &t& prises, Une 
bonne reconstruction de 1' image aux portions de 
recouvrement impose de combiner entre el les les mesures 
associ€es a chacun des details pour des vignettes 

20 diffgrentes, ce qui est impossible a faire directement 
puisque leurs hauteurs sont inconnues en g6nSral. Si on 
choisit par exemple d'associer les rayons distants de 
l'£cart 6 pour reconstruire les portions de 
recouvrement , les details 3. hauteur de 5a seront rendus 

25 correctement, mais les details presents k d'autres 
hauteurs ne pourront pas l'Stre. La combinaison des 
vignettes produira alors du flou et un grossissement 
inexact de ces autres details. 

Un procede de ce genre a cependant d§j& 6t§ 

30. propose- daiis 1'art. L' image, est reco'nstruite en 
choisissant une reconstruction exacte des portions de 
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recouvrement & des hauteurs deterrainees , Ik oil les 
details iraportants, et notamment les os pour une 
radiographie du corps, sont susceptibles de se trouver. 
Pour obtenir des bons r^sultats, il faut done connaitre 
5 d'abord la hauteur de ces details. Des coefficients de 
pond^ration peuvent privilggier les r^sultats de l'une 
ou 1' autre des vignettes selon la position consid£r6e 
sur la portion de recouvrement. La restitution des 
autres details de 1' image est sacrifice. 

10 Un autre proced£ connu consiste k calculer 

des correlations entre les portions de recouvrement des 
diff^rentes vignettes pour ^valuer 1'gcart (6, 7 ou 
autre) des rayons a associer pour faire la synthase des 
portions de recouvrement. Les correlations dependent de 

15 details pr6pond£rants presents sur les deux portions de 
recouvrement et provenant d'un merae endroit de l'objet 
2. La reconstruction de 1' image est accomplie hauteur 
de ces details preponderants et elle est bonne, si du 
moins ces details existent ; mais comme dans le proc€d€ 

20 precedent, les details situes aux autres hauteurs 
seront mal rendus. 

II faut aj outer que des conflits de hauteur 
peuvent apparaitre si les portions de recouvrement sont 
nombreuses, et en particulier avec un rayonnement 

25 conique ou les portions de recouvrement concernent tout 
le perim^tre des vignettes. Deux portions de 
recouvrement sur deux c6t£s d/une vignette pourront 
Stre reconstruites independamment k des hauteurs 
diff£rentes, tout en ayant une intersection pour 

30 laqiielle on sera erabarrasse a choisir une hauteur de 
reconstruction. 
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Un proc6d6 de reconstruction plus correcte 
d'une image radiographique est propose avec 
1' invention. II est fond6 sur une discretisation 
g^n^rale de 1' objet en volumes (voxels) definissant des 
5 hauteurs de reconstruction, et des combinaisons des 
valeurs d' attenuation estim£es sur chacun des volumes 
aux diff§rentes hauteurs de reconstruction pour 
ameiiorer 1' image globale, sans priviiegier 
n§cessairement une hauteur de reconstruction. 

10 De fagon plus detailiee, 1' invention 

concerne gene ralement un proc^de de reconstruction 
d'une image radiographique d'un objet traverse par un 
rayonnement divergent subissant une » attenuation, le 
rayonnement occupant des positions successives ayant 

15 des portions de chevauchement et 1 1 attenuation etant 
mesuree par un reseau de detect eurs, sur lesquels le 
rayonnement se projette et dormant des vignettes de 
1' image respectivement associees aux positions du 
rayonnement et comprenant aussi des portions de 

20 chevauchement, le procede comprenant une comb inai son de 
vignettes pour reconstruire 1' image, ainsi que les 
etapes suivantes : 

- discretiser l 1 objet en voxels definissant 
des hauteurs de reconstruction, 

25 - associer chaque voxel & au moins un 

detecteur respectif du reseau sur lequel le rayonnement 
se projette apres avoir traverse ledit voxel, 

attribuer une valeur d' attenuation a 
chaque voxel d»apr£s les valeurs* mesurees par ledit 

30 detecteur associe, 
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- et combiner les valeurs d 1 attenuation des 
voxels aux differentes hauteurs de reconstruction pour 
obtenir une image bidimensionnelle. 

Dans une des formes de 1' invention, la 
5 valeur d' attenuation attribute a chaque voxels est 
egale a la somrae des valeurs mesurees par ledit 
detecteur associe, divis^e par le nombre des vignettes 
que contribuent S. donner ledit detecteur- associe et par 
la longueur de chaque voxel qui a ete traversSe, et les 

10 valeurs d 1 attenuation des voxels sont combinSes par une 
cotnbinaison num6rique sur des groupes des voxels 
superposes aux differentes hauteurs de reconstruction. 
Et dans une autre de ses formes, la valeur 
d' attenuation attribute a chaque voxel est obtenue par 

15 retroprojection iterative des valeurs d' attenuation 
mesurees par les detecteurs, des valeurs provisoires 
etant attribuees aux voxels et corrigees aprds avoir 
ete projetees sur les detecteurs, en calculant des 
differences entre des sommes et les valeurs provisoires 

20 sur des lignes de projection aux valeurs mesurees par 
les detecteurs sur lesdites lignes de projection, et en 
repartissant les differences sur lesdites lignes de 
projection pour corriger les valeurs provisoires. On 
ef f ectuera aussi une combinaison numerique des . groupes 

25 de voxels. 

L' invention sera maintenant decrite 
compietement en liaison aux figures, dont la figure 1 
schematise le procede de realisation des vignettes, la 
figure 2 illustre le probl£me de reconstruction a une 
30 hauteur arbitraire, la figure 3 illustre les elements 
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explicatifs de 1' invention, et les f igures 4 et 5 sont 
des organigratnmes de deux modes du proced#. 

On passe au commentaire de la figure 3. 
L' ob jet 2 est discretise en volumes 61ementaires ou 
5 voxels 8 qui definissent des hauteurs de reconstruction 
11. lie rayonnement passe a travers les voxels 8 par des 
rayons 9i et 9 j , qui sont plusieurs pour les voxels 8 
appartenant aux portions de recouvrement, et qui sont 
originaires de positions respectives li et lj de la 

10 source 1 et se projettent sur des detecteurs respectifs 
lOi et lOj qui leur sont associes pour les positions 
correspondantes 3i et 3 j du r£seau de detecteurs 3 . Les 
detecteurs 10 mesurent des attenuations des rayons 9 a 
travers tout 1'objet 2, et done k travers tous les 

15 voxels tels que 8 qu'ils traversent. Dans la pratique, 
les voxels 8 se projettent sur une surface qui- peut 
englober plusieurs detecteurs 10 completement , et 
d'autres partiellement . Le systSme est calibre pour 
associer a chaque voxel 8 les detecteurs 10 sur 

20 lesquels il se projette et repartir sur eux les 
proportions de son attenuation. Nous n'evoquerons pas 
ici ces techniques de calibration, qui sont tout & fait 
usuelles dans 1'art, et considererons des voxels 8 se 
projetant completement sur un detecteur 10 unique le 

25 long d'un rayon de projection unique par souci de 
simplicite des explications. 

D'aprSs la figure 4, le procede commence 
par une etape A de discretisation generale de 1'objet 2 
en voxels 8 dont les couches definissent les hauteurs 

30 de reconstruction 11 de 1' image.* En pratique, les 
hauteurs de* * reconstruction 11 seront assez peu 
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notnbreuses et les voxels 8 plut6t parallelepip§diques, 
allonges en hauteur, que cubiques. lies etapes suivantes 
B et C consistent & se placer & une hauteur de 
reconstruction 11 et un voxel 8 determines . On 
5 recherche alors les rayons tels que 9i et 9j passant 
par le voxel 8 consid^re, et les ditecteurs 10 i et 10 j 
de projection desdits rayons sur le reseau de 
dStecteurs 3, a l'£tape D. L'etape suivante E consiste 
k lire la mesure d' attenuation des rayons 9i et 9j sur 

10 les d^tecteurs lOi et 10 j . A l'etape suivante F, une 
moyenne de ces attenuations est faite, du moins pour 
les voxels 8 appartenant aux portions de recouvrement 
et qui sont done traverses par au moins deux rayons 9, 
Une pond£ration peut Stre appliquee aux diffSrentes 

15 mesures, en accordant par exeraple plus de poids a 
eel les qui proviennent de rayons sensiblement 
verticaux, ce qui am^liore en particulier 1 1 image aux 
positions de de recouvrement des vignettes* Les Stapes- 
C a E ou F sont ensuite r£p£t£es pour tous les volumes 

20 de la couche consid£r£e ; apr&s quoi, a l'6tape G, une 
image de l'objet 2 est reconstruite. 

Cette image est une image de 1' ensemble de 
l'objet 2, et pas seulement une image de coupe a la 
hauteur considSr^e, puisque les attenuations mesurles 

25 par les dStecteurs 10 le long des rayons 9 orit £t£ 
suppos^es concentrees aux voxels 8 de la couche a cette 
hauteur. 

Ensuite, on revient a l'etape B pour 
reconstruire l'objet 2 a une autre hauteur, et le cycle 
30 des etapes C a G recommence avec les voxels 8 de la 
couche associee. 
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Quand les images de I'objet 2 ont £te 
reconstruites a toutes les hauteurs, elles sont 
combinees 3. I'etape H avec l'espoir d'obtenir une image 
plus exacte. Plusieurs procedes peuvent Stre envisages. 
5 Le plus simple consiste peut-gtre £ faire des moyennes 
des images sur des colonnes 12 (It la figure 3) de 
voxels 8 empiies appartenant a des couches diff ^rentes, 
avec eventuel lenient une ponderation pour favoriser les 
couches les plus representatives. Eventuel lement, on 

10 peut choisir une seule des images qu'on juge meilleure 
que les autres, ou un assemblage de plusieurs des 
images aux endroits qu' elles representent le mieux. 
Toutes ces methodes devraient donner de meilleurs 
resultats que celles de l'art anterieur qu'on a 

15 decrites auparavant . 

On ne fera que mentionner certains procedes 
de correction qui sont usuels dans l'art et qui ne sont 
pas affectes par 1' invention. 

Le rayonnement diffuse peut tout d'abord 

20 §tre soustrait des mesures avant d' exploiter celles-ci. 
Plusieurs procedes existent pour faire cette 
soustraction, dont le plus simple, qu'on donne a titre 
d'exeraple, est peut-Stre de proc€der & une mesure 
supplementaire oil un ecran est intercaie entre I'objet 

25 2 et le rgseau de detecteurs 3 en raasquant certains des 
detecteurs 10. Les detecteurs 10 masques ne sont pas 
touches par le rayonnement direct des rayons 9, mais 
seulement par le rayonnement diffuse, qui est alors 
mesure par ces detecteurs et qui peut etre deduit par 

30 des interpolations pour les autres detecteurs. 
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Les attenuations d'un rayonnement peuvent 
en g6neral etre exprimees par un coefficient 
multiplicateur du rayonnement initial Io inf^rieur a 
1' unite" et egal a oil 1 est la longueur 

5 d' attenuation et \i le coefficient d' attenuation 
caract^ristique du materiau, et qui est en general la 
valeur qu'on recherche pour reconstruire l f image. Les 
de"tecteurs 10 mesurent directement le rayonnement I qui 
n'a pas £te absorbe par l'objet 2 et qui est egal a 

10 I 0 e"* 1 ; on peut en deduire le produit fXl, puis la 
valeur de \l si on divise les valeurs de fll par les 
longueurs de travers6e de l'objet 2 par les rayons 9, 
apr§s les avoir estimees par une autre mesure ou les 
avoir evaluees geomStriquement . 

15 Un autre mode de realisation de 1' invention 

sera maintenant dScrit au moyen de la figure 5. Aprds 
une etape de discretisation J semblable a celle A de la 
realisation pr6c6dente, on effectue au mieux une 
division en blocs & l'Stape K. En effet, la resolution 

20 qui va Stre entreprise peut devenir malaisee si le 
systdme consid^re est trop volumineux. En pratique 
chaque bloc pourra comprendre les voxels 8 associe*s S 
une vignette. Qu'une division en blocs soit faite ou 
non, le probleme £ r^soudre peut s'exprimer par p = Mx 

25 oil x dSsigne les inconnues, c'est-£-dire les 
attenuations aux voxels 8, p d^signe les projections de 
ces valeurs, c'est-^-dire les raesures par les 
de*tecteurs 10, et enfin M designe la matrice de 
projection. Les coefficients mij de la matrice M 

30 representent la contribution d'un voxel 8 d'indice j a 
la projection suivant le rayon 9 d'indice i, et peuvent 
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en general §tre approchEs par la longueur traversee par 
ce rayon dans ce volume. 

L'Etape suivante L est une evaluation de 
1' attenuation aux voxels 8 du bloc considere. La 
5 premidre evaluation peut dtre arbitraire, par exemple a 
des valeurs nulles. Pour chacun des voxels 8 on 
recherche, d'apr^s 1'etape M, le detecteur 10 qui lui 
est associE par le rayon 9 qui le traverse, 
1'etape D de la realisation prEcEdente. L'Etape 

10 suivante N est une lecture des raesures des dEtecteurs 
10 semblable a 1'etape E . La determination des rayons 9 
de projection permet de procEder 3. une Evaluation des 
valeurs projetEes d' attenuation £ l'Etape O, c'est-t- 
dire qu'on procSde au calcul Mx pour Evaluer p. En 

15 soustrayant ces valeurs EvaluEes des projections aux 
valeurs reelles, mesurees, des mEmes projections, on 
determine I'erreur commise dans 1' evaluation des 
valeurs projetees a 1'etape P. 

L'Etape suivante Q est une rEtroprojection 

20 de cette erreur dans les voxels 8 de l'objet 2 afin de 
corriger les valeurs EvaluEes de 1' attenuation. 
ConcrEtement, on procede en executant la formule 

et x (q> sont des Evaluations successives de 1' attenuation 
25 aux voxels 8 du bloc ; A, (q) est tan coefficient de 
relaxation permettant- de ne pas se diriger trop vite 
vers une solution qui ne correspond qu'aux premiers 
blocs et qui est corapris entre 0 et 2 ; ce coefficient 
n'est d'aille.urs pas uniforme dans les blocs maj.s peut 
30 avantageusement Etre plus ElevE pour les rayoris 
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sensiblement verticaux, ou perpendiculaires aux 
detecteurs 10, afin de leur dormer une importance de 
pondSration plus grande, comme dans la realisation 
pr^cgdente ^loc est la transpos^e de la matrice M 
5 pour le bloc considSre ; le terme au d§nominateur est 
un terme de normalisation ; enfin, les termes entre 
parentheses repr^sentent l'erreur calculee a l'etape P. 

On procede de mSrae pour le bloc suivant , en 
recommengant le cycle de l'etape K & l'etape Q, puis on 

10 revient au premier bloc pour une nouvelle iteration, 
jusqu'a ce que les attenuations evaluees aient converge 
vers une solution, ce qu'on exprime par l'etape R. Les 
voxels 8 inclus dans les portions de recouvrement de 
I'objet 2 ont ete traites de la meme fagon que les 

15 autres, en subissant simplement des iterations plus 
norabreuses, pour chacun des blocs auxquels ils 
appartiennent . 

On dispose alors d'une image 

tridimensionnelle de l'objet 2 ; une image 

20 bidimensionnelle de bonne qualite peut etre obtenue par 
vine combinaison des valeurs obtenues, qui consiste £ 
combiner les valeurs d' attenuation sur les colonnes 12 
de volumes 8 empiies. 

Le proced§ de 1' invention permet de 

25 concilier une bonne qualite de restitution des details 
importants de l'objet etudie avec une bonne qualite 
d' ensemble de 1' image. II est possible d'obtenir des 
images dont la resolution est analogue au pas des 
detecteurs 10 . 

30 On s'est place dans la situation usuelle oft 

le reseau 3 de detecteurs accompagne le mouvement du 
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rayonnement 4, mais le proc^de pourrait §tre appligu6 
sans changement avec un rSseau de dStecteurs inunobile 
sous l'objet 2 et dont la surface s'§tendrait a toutes 
les vignettes de projection. 



10 



WO 2004/059578 



PCT/FR2003/050195 



14 



REVEND X CAT I ONS 
1) Proc§d6 de reconstruction d'une image 
radiographique d'un objet traverse par un rayonnement 
divergent subissant line attenuation/ le rayonnement 
5 occupant des positions successives (4) ayant des 
portions de chevauchement et 1 1 attenuation etant 
mesur^e par un r£seau (3) de detecteurs (10) , sur 
lesquels le rayonnement se projette et donnant des 
vignettes de I 1 image respectivement associees aux 
10 positions du rayonnement et comprenant aussi des 
portions de chevauchement , le procede comprenant une 
combinaison de vignettes pour reconstruire 1" image, 
ainsi que les Stapes suivantes : 

discretiser l 1 objet en voxels (8) 
15 definissant des hauteurs de reconstruction (11), 

associer les voxels & au moins un 
detecteur respectif du r^seau sur lequel le rayonnement 
se projette apr£s avoir traverse ledit volume, 

attribuer une valeur d 1 attenuation & 
20 chaque voxel d'apres les valeurs mesurees par ledit 
detecteur associe, 

- et combiner les valeurs d 1 attenuation des 
voxels aux diff£rentes hauteurs de reconstruction pour 
obtenir une image bidimensionnelle. 
25 2) Procede de reconstruction d'une image 

radiographique selon la revendication 1, caracterise en 
ce que la valeur d' attenuation attribuee a chaque 
volume est egale k la somme des valeurs mesurees par 
ledit detecteur associe, divisee par le nombre des 
30 vignettes que- contribuent a donner ledit detecteur 
associe, et les valeurs d' attenuation des voxels sont 
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combinees par une cornbinaison num§rique sur des groupes 
(12) des voxels superposes aux diff-erentes hauteurs de 
reconstruction . 

3) Procede de reconstruction d'une image 
5 radiographique selon la revendication 1, caracterise en 

ce que la valeur d 1 attenuation attribute a chaque voxel 
est obtenue par retroprojection iterative des valeurs 
d 1 attenuation mesurles par les detecteurs (10), des 
valeurs provisoires etant attributes aux voxels et 

10 corrigees apres avoir §te projetees sur les detecteurs, 
en calculant des differences entre des sotnmes des 
valeurs provisoires sur des lignes de projection et les 
valeurs mesurees par les detecteurs sur lesdites lignes 
de projection, et en retroproj etant les differences sur 

15 lesdites lignes de projection pour corriger les valeurs 
provisoires. 

4) Procede de reconstruction d'une image 
radiographique selon la revendication 2 ou 3, 
caracterise en ce que les valeurs d f attenuation des 

20 volumes sont combinees numeriquement sur des groupes 
(12) des volumes superposes aux differentes hauteurs de 
reconstruction. 

5) Procede de reconstruction d f une image 
radiographique selon I'une quelconque des 

25 revendications precedentes, caracterise en ce qu'il est 
applique & une osteodensitometrie . 
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